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While the accelerated expansion of the Universe is by now well established, an underlying scalar 
field potential possibly responsible for this acceleration remains unconstrained. We present an at- 
tempt to reconstruct this potential using recent SN data, under the assumption that the acceleration 
is driven by a single scalar field. Current approaches to such reconstructions are based upon sim- 
ple parametric descriptions of either the luminosity distance or the dark energy equation of state 
(EOS). We find that these various approximations lead to a range of derived evolutionary histories 
of the dark energy equation of state (although there is considerable overlap between the different 
potential shapes allowed by the data). Instead of these indirect reconstruction schemes, we dis- 
cuss a technique to determine the potential directly from the data by expressing it in terms of a 
binned scalar field. We apply this technique to a recent SN dataset, and compare the results with 
model-dependent approaches. In a similar fashion to direct estimates of the dark energy equation 
of state, we advocate direct reconstruction of the scalar field potential as a way to minimize prior 
assumptions on the shape, and thus minimize the introduction of bias in the derived potential. 

PACS numbers: 98.80.Es, 97.60.Bw, 98.80.Cq 



I. INTRODUCTION 

Distance estimates to Type la supernovae (SNe) are 
currently a preferred probe of the expansion history of 
the Universe and have led to the discovery that the 
expansion is accelerating It is now believed that a 
mysterious dark energy component, with an energy den- 
sity ~70% of the total energy density of the universe, 
is responsible for the accelerated expansion Q. While 
the presence of acceleration is now well established by 
various cosmological probes, the underlying physics re- 
mains a complete mystery. As the precise nature of the 
dark energy has profound implications, understanding its 
properties is one of the biggest challenges today. 

With the advent of large surveys for Type la super- 
novae, such as the Supernova Legacy Survey (SNLS) f26j 
and Essence [13], among others, it is hoped that we will 
study details of the expansion, and thereby elucidate the 
physics responsible for the acceleration. Under the as- 
sumption that the dark energy is due to a single scalar 
field rolling down a potential, several studies have con- 
sidered how future data might be used to reconstruct the 
potential, either based on various analytical descriptions 
of the luminosity distance Q , or through specific assump- 
tions about the potential, such as a polynomial function 
in the scalar field It is already well established that 
certain parametric descriptions of the distance lead to 
biased estimates for the dark energy cquation-of-state 
(EOS) and the potential Q. While improved paramet- 
ric forms of fitting functions have been suggested [1, [l^l , 
it is unclear how to select an optimal approach for re- 
constructing the dark energy scalar field potential from 
SN distances (for a review of various possibilities, see 
Ref. [iH). 



In this paper we discuss issues related to potential and 
dark energy EOS reconstruction by making use of a re- 
cent set of SN data from the SNLS survey [l^l . The sam- 
ple includes 73 high redshift SNe complemented with a 
sample of 44 nearby supernovae [l2|. We compare and 
contrast a variety of methods to reconstruct the potential 
and the dark energy EOS. We write the luminosity dis- 
tance either as a simple polynomial expansion in redshift, 
or as a Fade approximation [l^ (which avoids some of the 
known problems in the polynomial expansion when tak- 
ing derivatives 0,11, El)- In addition to approximating 
the luminosity distance, we also explore two approxima- 
tions to the EOS: w{z) = wq + Wa{l — a) [H, [l^ and 
w{z) = Wo - aln(l + z) [1]. 

Based on our model reconstruction of the potential, 
we find that while there is significant overlap of the al- 
lowed V{(j)) region favored by each of the four reconstruc- 
tion methods, the models give rise to different histories 
for the EOS, especially within the two parameter plane, 
w-w' (the EOS parameter, w, and its time derivative, 
w' = dw/d\na, as functions of redshift We argue 

that existing parametric fitting functions for either dis- 
tance or the EOS lead to biased reconstructions of the 
potential. In the literature, however, there exist model- 
independent ap pro aches to the reconstruction of the dark 
energy density [Tsj and the EOS which bin the pa- 
rameters directly as a function of redshift, with the num- 
ber and width of the bins determined by the statistical 
quality of data. These estimates can also be arranged 
to be uncorrelated [Toj , allowing unique insights into the 
evolution without being subject to prior assumed red- 
shift dependencies. Here we suggest a similar model- 
independent approach to the reconstruction of the scalar 
potential from SN data. Instead of utilizing a polynomial 



2 



expansion for the potential which assumes a Hmited 
range of models (once the expansion is truncated at a 
certain order), we propose a binning scheme for the po- 
tential that can be applied to data with a minimal, and 
easily controlled and understood, number of assumptions 
for the potential shape. 

The paper is organized as follows: In the next Section 
we review techniques for reconstructing the scalar-field 
potential from SN distances. We also reconstruct the 
EOS as a function of the redshift, and use this to study 
the w-w' plane (which has been advocated as a way to 
characterize the underlying potential responsible for the 
dark energy component by separating the regime into 



"freezing" and "thawing" potentials [17|; see, also |20| 
for a Monte-Carlo exploration). In Section III we explore 
the impact of different paramctcrizations on the derived 
evolutionary histories. While we observe these differences 
with ~ 115 SN data points, future large SN datasets may 
lead to apparently inconsistent results. In Section IV, fol- 
lowing the approach to model-free estimates of the dark 
energy EOS [l9[, we present a model-independent esti- 
mate of the scalar field potential. Wc conclude with a 
summary of our main results in Section V. 



II. POTENTIAL VIA PARAMETRIC FORMS 

For this study we make use of SN data from SNLS [l^l . 
Due to complications related to independent data sets 
(e.g., differing calibration, color correction, extinction 
correction, etc.), we do not attempt to increase the sam- 
ple size by combining other SN datasets. The measure- 
ments from Ref. [13 present the quantity fiB = tub — M 
for 117 SNe, with 73 of these at redshifts greater than 0.2 
[28t This distance modulus is related to the luminosity 
distance through /is = Slogj^gdi, while the luminosity 
distance is related to the comoving radial distance via 
dL = c(l + z)r{z)/HQ, where r(z) = dz'/H{z') with 
H{z) the expansion rate of the Universe. When model 
fitting the data, we fix M = 19.3 ± 0.03 to the value 
determined by SNLS. We take the central value; further 
uncertainty will be incorporated into Cint, as discussed 
below. 

In our reconstruction of the potential, we describe r{z) 
through two parametric forms widely used in the litera- 
ture. First, we expand r{z) as a simple power-law Q 
such that 



r{z) = z + a2z'^ + a^z^ + a^z^ . 



(1) 



Note that the coefficient of the first order term is ex- 
actly one. Since this polynomial expansion has known 
problems when estimating the derivatives of r{z) (e.g.. 
Figure 3 of Ref. and also Ref. [§]), we also consider a 
Fade form for r(z) with Ref. [l3l |: 
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FIG. 1: Hubble diagram for 115 Type la SNe used in the 
present analysis. The error bars are on /is only. We also 
include an additional constant error, aint = 0.13, to account 
for the SN intrinsic dispersion. For reference we also plot 300 
curves drawn uniformly from the 2cr consistent likelihood fits 
to the data using the Taylor expansion with r(z) = z + a2Z^ + 



such that as 2: ^ 0, r[z) — s- z. In this form, using r{z — > 
00), one can additionally constrain the parameters with: 
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In addition to the two fitting forms for r(z), we also de- 
termine r(z) through model parameterizations for w{z), 
including w(z) = wq + {1 — ajWa [H, [3 and w{z) = 
Wo + aln(l -I- z) [1]. Since from w it is possible to deter- 
mine the distance, these approximations allow us to once 
again reconstruct the dark energy potential. 

In each of two parametric descriptions of r{z) we have 
three free parameters. We parameterize wlz) with two 
parameters, and include Qm as a third free parameter 
(under the assumption of a flat universe; weakening this 
assumption significantly degrades our ability to measure 
anything about the potential with existing data). When 
showing results related to potentials or EOS as a function 
of redshift, we take a prior on f2,„ such that the prob- 
ability is Gaussian with a mean of 0.25 and a standard 
deviation given by cr = 0.05 [3. In each case, to obtain 
the join likelihood distribution of the parameters given 
the data, we perform a likelihood analysis: 



N 



(4) 



where, following Ref. [1^, in addition to statistical un- 
certainty in fiB we include an additional Gaussian uncer- 
tainty, (Tint — 0.13, representing the intrinsic dispersion 
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of SN absolute magnitudes. M . We ignore complications 
related to covariances in the Hubble diagram, either due 
to effects related to calibration [2^ or fundamental lim- 
itations such as gravitational lensing correlation of SN 
flux [13] or peculiar velocities [2^. The posterior prob- 
ability distribution is taken to be P{pi\^) oc e~5X ivi) ^ 
and we marginalize the likelihood over the uncertainty in 
rim, assuming a Gaussian prior distribution. 

Once the joint probability distribution for parameters 
is determined, we sample the Icr and 2a range allowed 
by these parameters to draw a fixed (> 600) number of 
independent r{z) curves consistent with the data. For 
each of these distance curves, ri{z), we obtain the scalar- 
field potential, in dimcnsionless units such that V{4') = 
^(0)/Pcrit = F/(3i/2/8^G), through Ref. i 
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where f — H^r. For each of the ri{z) estimates, we also 
randomly draw f2m from a Gaussian prior distribution 
as described above. The mapping between z and (f)^ the 
scalar field value, is obtained through 



d4> dr I dz 

dz ~ (1 + z ) 

1 {l + z)d^f/dz^ 
An {df/dzY 
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FIG. 2: The normalized quintessence potential V{(f>) vs 
(5) 4>/mpi. The shaded region is allowed at the 2a confidence 
level when using the Taylor expansion for r{z). The red 
lines mark the same when using a Fade approximation to 
the distance. The orange and blue lines are for the cases 
where w{z) is parameterized by w{z) = -I- aln(l -f z) and 
w{z) = Wo -\- Wail — a), respectively. Here = corresponds 
to z = 0, while > 0.1 generally corresponds to z > 1 (de- 
pending on d(j>/dz). The points with error bars show the la 
(solid) and 2a (dashed) model-independent estimates of the 
potential described in Section IV (see equation (8)). While 
there is considerable overlap in the allowed region, there are 
also significant differences in terms of the redshift evolution 
of the EOS. 



where cf) = (p/m-px. Furthermore, for models where we 
parameterize r(z), wc can also extract the dark energy 
EOS as 



w{z) 



1 + z 3an(l + zf + 2{d^ri/dz^)/{dn/dzf 
~3 n„^{l + z)3 - (drjdz)-^ 
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When selecting models associated with scalar fields, we 
require that d<j)/dz > 0, such that w > —1. Even in the 
case of w{z) parameterizations where model fits allow 
w < — 1, we ignore w{z) below this value as single scalar- 
field models do not naturally give rise to such EOS. 



III. BIASES IN MODEL-DEPENDENT 
ESTIMATES 

In Fig. 1 wc show the Hubble diagram for the 115 data 
points from Ref. (l^ used in this analysis. For reference, 
we also plot ^ 300 distance curves which arc 2cr consis- 
tent curves drawn from the likelihood distribution for pa- 
rameters under the Taylor expansion for r{z). The best- 
fit model with this parameterizations has a chi-square 
value of 113.1 with 112 degrees of freedom. Note that 
in Ref. [l2l, tTjnt is tuned so that = 1 for the best- 
fit model under standard-ACDM cosmological fits to the 



data. We use their best-fit value, tJint = 0.13, and do not 
take this intrinsic uncertainty as an additional free pa- 
rameter in our modelling. The exact value of the intrinsic 
dispersion does not impact our comparison of different 
approaches to the reconstruction of the quintessence po- 
tential. It is to be emphasized that all of our parameter- 
izations of either distance or the EOS yield comparable 
values for the best-fit model. This suggests that all 
four of the reconstruction methods outlined above are 
indistinguishable within the redshift range considered. 

As discussed in the previous section, for each of the 
four parameterizations we determine a best-fit r(z) to the 
SN data. We then Monte-Carlo generate models within 
2a of this best-fit, generating over 600 instances of w{z) 
and V{(j}), all of which are consistent with the underlying 
SN dataset at the 2a level. In Fig. 2 we show the poten- 
tials reconstructed from each of the four methods, with 
the bands encapsulating 95% of the distribution of the 
individual models. Due to the behavior of the Taylor ex- 
pansion at high 2, and the fact that we do not restrict the 
coefficients of the polynomial expansion to follow a flat 
universe, this parametrization gives rise to a large range 
of acceptable potentials which satisfy the data. The Fade 
parametrization of r(z), as well as the w(z) models, sig- 
nificantly improve the constraints on allowed potential 
shapes. This is because the parameters in the Fade ap- 
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FIG. 3: The dark energy EOS, w{z), as a function of red- 
shift. The curves show the 2a allowed values and correspond 
to the potentials shown in Fig. 2. Note that if we impose 
d(j)/dz > 0, then w > —1, as single scalar-field models do not 
lead to ui < —1. However, as shown, most of the parameter- 
izations allow the region below ui < — 1. When constructing 
the potentials shown in Fig. 2, we apply the condition that 
d(j>/dz > 0. Note that the w{z) parameterizations, with two 
free parameters, are the most restrictive parametrization in 
the regime z < 0.3. Over the redshift range probed, the dif- 
ferent parameterizations generally agree with each other. The 
plotted error bars show the la and 2a errors of 'Wi{z) when 
the EOS is subdivided into three bins in redshift, with 
directly measured from data and no restrictions on its values. 
A Gaussian prior has been taken on Sim with a one sigma 
uncertainty of 0.05 with w{z) parameterizations. 



proximation are additionally constrained to satisfy crite- 
ria related to the behavior of r{z) as z — *■ oo, as well as 
by the assumption of a flat universe When fitting 

the w{z) parameterizations to the data, we were able to 
impose a prior on Qm based on existing cosmological in- 
formation (this was not possible when using r{z) fitting 
functions) . While we find some overlap in the 2cr allowed 
region in the (f>-V(4>) plane between the four approaches, 
there are also noticeable inconsistencies. Analysis of an 
identical dataset with different approximations to V or 
w lead to differing resulting best-fit potentials. 

The differences related to potential shapes between the 
four methods are best captured in terms of evolutionary 
histories for the dark-energy EOS. In Fig. 3 we summa- 
rize the best-fit w{z) results for each of the four recon- 
struction techniques. Note that some of our parameteri- 
zations allow w{z) < — 1, but due to our assumption that 
the dark energy arises from a scalar-field potential where 
w{z) is always expected to be greater than -1, we re- 
strict the allowed parameter space to be the region where 
w{z) > —1. Similarly to Fig. 2, we find considerable 
overlap between different reconstruction schemes in the 
w{z) versus redshift plane, with most models indicating 
that as the redshift is decreased, w{z) tends to values 



between -0.8 and -1.0 at z = 0. In terms of our direct 
w{z) parameterizations, with w{z) ^ wq + {1 — a)'Wa we 
find Wo = -1.12 ± 0.14 and Wa ^ 0.38 ± 0.49 at the 68% 
confidence level. In the case of w(z) = wq + aln(l -I- z) 
we find Wo = -1.08±0.11 and a = 0.35±0.75. As shown 
in prior studies [l3 |. w{z) parameterizations allow for a 
minimum w{z) region at a certain pivot redshift. For 
the dataset used here, this pivot redshift is at z ~ 0.12, 
and at the 2a confidence level we find that the pivot 
point satisfies —1.23 < Wp = w{z = 0.12) < —0.74, using 
w{z) = Wo + (1 — a)wa. It is important to note that all 
the parameterizations are consistent with a cosmological 
constant. 

In Figs. 2 and 3 we show the 2cr bands of best-fit mod- 
els to the data, under different parameterizations of the 
distance or dark energy EOS. In addition to this outer 
envelope, we are also interested in the distribution of 
the individual w{z) models within the 2a bands. We 
thus study the behavior of the models in the w-w' plane, 
which has been sug gest ed as a natural venue in which to 
distinguish models [ol. We Monte-Carlo 600 scalar po- 
tentials, l^((/>), and evolution histories, w{z), within the 
2a regime of the best-fit parameters for each of the four 
fitting functions. In Fig. 4 we plot w and w' sX z = 0.1 
and z = 0.5 for each Monte-Carlo model, with the scat- 
ter of points being 2a consistent with our underlying SN 
data set. 

Based on the evolutionary behavior of simple scalar- 
field models in the w~w' plane, it has been suggested that 
one can separate potentials into "thawing" and "freez- 
ing" regions, based upon their shapes [l3|. These regions 
are delineated in Fig. 4, for comparison with our individ- 
ual Monte-Carlo models. It is apparent that the different 
parametrization approaches yield separate, though often 
overlapping, regions within the w-w' plane. In addition, 
the models are not necessarily well-contained within the 
thawing or freezing regions, with a freezing model in one 
parametrization ending up as a thawing model in an- 
other, or with models ending up in between thawing or 
freezing, or well outside of either regime. Using generic 
numerical models for the potential shape, this behavior 
has also recently been highlighted in Ref . [20| . By apply- 
ing additional constraints on allowable potentials (espe- 
cially at high z), Ref. [l3l find much tighter confinement 
in the w~w' plane. 

Any statement regarding the shape of the scalar po- 
tential, as determined from data, is thus crucially depen- 
dent upon the underlying parameterizations. For exam- 
ple, for the Taylor expansion approach w' is largely neg- 
ative at z — 0.5, while it is positive at z = 0.1. Under the 
Fade approximation, w' is negative at both z = 0.1 and 
z = 0.5, with w tightly clustered (-0.9 < w < -0.8) 
at z — 0.1 and relatively unconstrained at z = 0.5. 
Although the reconstructed potentials show significant 
overlap (see Fig. 2), the distributions in the w~w' plane 
are less consistent among different parameterizations. 
Thus, while there is motivation from theoretical argu- 
ments for using the w-w' plane for potential recognition. 
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FIG. 4: w versus dw/dlna. The data points show the EOS 
and its time derivative for 600 model potentials uniformly 
drawn at the 2a confidence level at redshift of 0.1 (filled sym- 
bols) and 0.5 (open symbols). The cyan, red, orange, and 
blue data points show potentials selected under the Taylor 
expansion. Fade approximation, wi^z) = uiq + a ln(l + z) and 
w(yZ) — wo+Wail — a) model fits to the data, respectively. For 
comparison, we also plot the thawing (dashed) and freezing 
(solid) potential regions, following Ref. [ij- There are con- 
siderable differences in w and dw/d\na values, and the evolu- 
tion captured by two redshifts, between the four approaches. 
These values do not satisfy the expectations under simple 
model criteria for the dark energy potentials, though all of 
these potentials, and the w{z) curves, are consistent with the 
data at the 2a confidence level. 

there is no obvious, parametrization independent way to 
convert distance data to constraints in this plane. 

The differences seen in Fig. 4 are attributable to the 
difFerent parametric forms used to approximate the dis- 
tance or the dark energy EOS. To paraphrase our results: 
you get out what you put in. Furthermore, the fitting 
forms to both the distance and the EOS are motivated 
by their ability to fit data, and possess no clear physi- 
cal motivation. Fig. 4 thus emphasizes the need for an 
approach which makes minimal assumptions about the 
underlying potential, thereby maximizing the measure- 
ment of a completely unknown scalar field. The less we 
assume about the potential, the more powerful the en- 
suing measurement of its shape. Such an approach is 
presented in the following section. 

IV. MODEL-FREE ESTIMATES 

Thus far we have discussed results based on assumed 
parameterizations for either distance or dark energy 
EOS. These parameterizations lead to conclusions that 
are subject to the assumed parameterizations. It is de- 
sirable to make model independent estimates of dark en- 
ergy. In the case of the EOS w(z), one could approach 
this by binning w{z) in redshift [19| . Applying this to our 
SNLS dataset, we evaluate w{z) over three bins in red- 



shift, < zi < 0.25, 0.25 < Z2 < 0.6, and 0.6 < < 1.0, 
assuming w{zi) constant in each bin. 

The resulting best-fit to the SN data is shown by the 
data points with Icr and 2cr error bars in Fig. 3. We find 
w{zi) = -0.88 ± 0.28 and w{z2) = -1.02^^^^ with no 
useful constraint for w{z) in the Z3 bin. Although these 
bins arc correlated at the 10% level, it is possible to decor- 
relate the binned Wi (z) estimates following the approach 
of Ref. [1^. While only three bins are attempted here, as 
SN sample sizes increase, one can consider larger num- 
bers of bins, each narrower in redshift. The estimates 
shown in Fig. 3 are consistent with estimates based on 
both fitting functions to the EOS, w{z) = wo + Wa(l — a) 
and w{z) ~ wo + a\n{l + z). As discussed in Ref. [101, the 
binned estimates capture the dark energy EOS with min- 
imal prior assumptions on the parameterization. This is 
expected to maximize the information one can extract 
from the data, while minimizing the introduction of bi- 
ases. 

As discussed and noted elsewhere 0, Q , the scalar- field 
potential reconstruction is also subject to prior assump- 
tions on the fitting form. To avoid biases and to make 
statements that arc not subject to assumed parametriza- 
tion, it is useful to directly construct the potential from 
data. Recent approaches in the literature consider fitting 
distance data to a potential expanded as a polynomial in 
the scalar field with y(0) = E,=o ^^'Z'' SMS- Since 
we are forced to truncate the expansion at low order (for 
example, at cubic order with existing data 0), the poten- 
tial is no longer arbitrary, but rather has a very limited 
range of possible shapes. 

Instead of assuming a specific family of shapes for 
V{(f>), we propose a "model- free" extraction of the po- 
tential directly from the data. We make two assump- 
tions about the scalar field potential: (1) that it is a 
piecewise continuous function, and (2) that its structure 
is "uniform" in the (j) range explored by the data. For 
(A^ — 1)A0 < (j) < NA(f>, we describe the potential as a 
function of the field with constant gradients, dV/d(j), over 
binned intervals, A0: 

N-l 

V{(l)) = Vo+Y,idV/dcl)),Aq)+i<l)-iN-l)Aq)){dV/d(j))N ■ 

(8) 

Assumption (1) above ensures continuity of V{(j)), which 
is necessary since one evolves the potential through the 
dynamic equation for the field as (f> + 3H(j) + dV/ dcj) = 
0; and discontinuities would lead to infinite deriva- 
tives. This requirement is unnecessary when consider- 
ing parameter-free estimates of the dark energy EOS, 
which is allowed discontinuous jumps in redshift. Both 
the constant-value and the constant-slope approaches to 
parameterizing the dark energy EOS lead to similar con- 
clusions [1^. Assumption (2) states that our bins in (j) 
are fixed width: Atp is a constant, independent of 4>. 
This assumption could be relaxed (e.g., finer bins near 
= 0), but this would lead to additional parameters, 
in addition to introducing model-dependent assumptions 
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into the analysis. The expansion of V{4>) in Eq. [8] ap- 
pears to make the least offensive assumptions possible, 
and thereby offers the basis with which to maximally 
constrain the full range of possible underlying potentials. 

We apply the above potential description to SNLS 
data following the same approach as Ref. Q , with three 
free parameters: Vq, {dV/d(j))i ior < (j> < 0.03, and 
{dV/d(t))2 for (f) > 0.03. The sizes of the bins are cho- 
sen by the range of (f> we are able to constrain, which 
is in turn related to both the redshift range of the SN 
dataset and the shape of the potential. Note that we 
take (j) = to coincide with z = 0. Instead of {dV/d(p)2, 
we convert the gradient to a data point at </> = 0.06, al- 
though we find only an upper limit, as this gradient is not 
strongly constrained by existing data. In Fig. 2 we show 
the estimated potential and error bars at the Icr and 2a 
level. The potential values allowed by the data are gener- 
ally consistent with other indirect reconstructions based 
on fitting forms for the distance or the EOS. While fit- 
ting forms lead to largely positive F(^) at </> > 0.05, our 
binned approach finds only an upper limit in this range. 

While we have described the potential with only three 
parameters, this can be straightforwardly generalized to 
additional bins as the statistics and quality of the SN 
samples improve. In addition, we make a minimal num- 
ber of assumptions regarding the potential, and thus are 
not biased for or against any particular shapes for the 
scalar field potential. The proposed approach is similar 
to the case where the EOS is binned and directly mea- 
sured from the data without specifying a model for the 
evolution. As SN data samples increase in size, we be- 
lieve such a model independent approach will become a 
powerful tool in extracting information about underlying 
scalar field potentials. 

V. SUMMARY 

We have presented a reconstruction of a single scalar- 
field potential using recent SN data from the SNLS sur- 
vey [1^. We have shown that reconstructions based on 
various approximations to the distance and the EOS lead 
to differing evolution histories of the dark energy EOS, 
particularly when the models are examined in the w-w' 
plane. In this plane the same data can lead to large 



movements in best-fit models, depending on the specific 
approximations to distance or EOS which are being uti- 
lized. Thus the underlying model assumptions lead to 
biases, compromising our ability to distinguish evolution- 
ary behaviors of the dark energy. At present the mod- 
els are only weakly constrained by the data, and thus 
this model-dependence, although apparent, is not a crit- 
ical failure. As the data improves, however, a model- 
independent approach will be essential to determining 
an otherwise completely unknown scalar-field potential. 

As an alternative to existing indirect reconstruction 
schemes, we have thus proposed a technique which es- 
tablishes the potential directly from the data, with only 
minimal assumptions about the underlying shape of the 
potential. We take the potential to be a binned scalar 
field, piecewise linear and continuous, but otherwise com- 
pletely arbitrary. Given the simplicity of these assump- 
tions, this potential is unlikely to introduce biases in the 
determination of a completely unconstrained, underly- 
ing potential. We have demonstrated this approach with 
current SN data, comparing the results to parameterized 
analysis. The ensuing constraints, although weaker, are 
expected to be robust and unbiased. It has been found 
that direct binning approaches to the dark energy EOS 
hold great promise for establishing model-independent 
measurements [l9| . We propose a similar approach to re- 
constructing the underlying dark energy scalar field po- 
tential, allowing us to make assumption-free statements 
about the nature of the completely unknown and mys- 
terious field potentially responsible for the accelerating 
expansion of the Universe. 
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